Domain wall dynamics in a two-component Bose-Mott insulator 
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We model the dynamics of two species of bosonic atoms trapped in an optical lattice within the 
Mott regime by mapping the system onto a spin model. A field gradient breaks the cloud into 
two domains. We study how the domain wall evolves under adiabatic and diabatic changes of this 
gradient. We determine the timescales for adiabaticity, and study how temperature evolves for 
slow ramps. We show that after large, sudden changes of the field gradient, the system does not 
equilibrate on typical experimental timescales. We find interesting spin dynamics even when the 
initial temperature is large compared to the super-exchange energy. We discuss the implication of 
our results for experiments wishing to use such a two-component system for thermometry, or as part 
of a cooling scheme. 



I. INTRODUCTION 



A major frontier in ultracold-atom physics is the search 
for novel low temperature (T) phases of matter. One well 
established success is the observation of Mott insulating 
shells of bosons in an optical lattice[l|, [2(- More recently 
Jordens et al. Q successfully observed the related Mott- 
metal crossover in a gas of fermionic 40 K atoms in an 
optical lattice. Much interest is now focussed on the mag- 
netic properties of these insulators, which set in at the 
super-exchange scale T ~ t 2 /U, where t is the tunneling 
amplitude, and U is the on-site interaction. For example, 
Trotzky et al. [J] observed physics of super-exchange in 
a double well setup. This physics has not yet been seen 
in an extended system. The major impediments to these 
studies center around cooling, thermometry and equili- 
bration. By considering an experimentally relevant setup 
involving two-component bosons, we are able to better 
characterize these impediments and explore ways to over- 
come them. We identify time-scales for adiabaticity, and 
the limits on the use of adiabatic demagnetization for 
cooling. 

The problem of cooling atoms in an optical lattice is 
particularly rich. It has been the focus of a large num- 
ber of theoretical works: for example M. Popp et at, B. 
Capogrosso-Sansone et al., J-S. Bernier et al. and T. 
L. Ho et al. 043 considered separating the system into 
entropy rich and poor regions, and subsequently evapo- 
rating away atoms in the entropy-rich regions. The lat- 
ter authors [9( also proposed a scheme where one reduces 
the entropy per particle in the region of interest by intro- 
ducing a defornration in the confining trapping potential 
as well as pointing out possible issues reg arding heat- 
ing while loading into an optical lattice [lfl|. Bezett and 
Blakie [ll[ have performed a study of the adiabatic path- 
ways for cooling fcrmions loaded into a three-dimensional 
optical lattice. J. Catani et al. 12[ proposed using 



species-selective trapping potentials to transfer entropy 
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from the target species to the auxiliary species. Addi- 
tionally, F. Werner et al. [l3[ have proposed using an 
analog of Pomcranchuk cooling in He-3 to adiabatically 
cool while loading a two-component Fermi gas in an op- 
tical lattice. 

Our study has relevance to all of these protocols. The 
key feature of each of these ideas is that they rely on adi- 
abaticity. Among the questions we address is: how slowly 
an external tunable parameter should be changed in or- 
der to maintain thermal equilibrium throughout? Not 
surprisingly, we find in the regime where super-exchange 
is important, the time-scale for adiabaticity scales as the 
inverse of the super-exchange energy. 

In addition to the issues of cooling and adiabaticity, 
a large amount of work has been done on thermometry 
[lj-[la]. Recently D. Weld et al. have suggested that a 
two-component Bose gas in the Mott regime can be used 
as an accurate thermometer. By displacing the poten- 
tials for the two species, one creates a domain wall. The 
width of that domain wall (w) is proportional to the tem- 
perature. It is this two-component Bose system that we 
study here. The question of thermometry is intimately 
connected to the question of adiabaticity and equilibra- 
tion. The system needs to be in thermal equilibrium for 
temperature to be well defined. 

Finally, our studies of the dynamics of this system fit 
in with ong oing works on ultra-cold atoms far from equi- 
librium J18H301 ]. Cold atoms systems are readily isolated 
from their environment, and their parameters may be 
rapidly tuned. This is largely unexplored territory, and 
organizing principles have only been found for a small 
subset of the possible phenomena. The response of this 
bosonic system to a changing field gradient is one of the 
more controlled avenues for exploring non-equilibrium 
physics. 

This paper is organized as follows: In Sec. II we dis- 
cuss the experiment of Weld et al. and introduce our 
model. In Sec. Ill we discuss the thermodynamics of the 
nrodcl and describe our numerical procedure for calculat- 
ing dynamics. In Sec. IV we determine the time-scales 
required for maintaining global equilibrium, and in Sec. 
V we discuss the idea of adiabatic demagnetization ap- 
plied to our existing setup. In Sec. VI we show that ef- 



fects of superexchange can be manifested even at larger 
temperatures. 



II. SPIN-GRADIENT THERMOMETRY 

In [3, Weld, Medley, Miyake, Hucul, Pritchard and 
Ketterle created a two-component Mott insulator out of 
two hyperfine states of 87 Rb in a cubic lattice. They 
worked well within the Mott regime using a lattice with a 
depth Vq = U.5E R (t/U = 0.024). At unity filling factor, 
this should be contrasted with lattice depth of 13. \Er 
(t/U = 0.034) for the supcrfluid to Mott-insulator tran- 
sition in three dimensions [17J. Here Er — -^-y is the 
lattice recoil energy where the lattice spacing is A, the 
atomic mass is m and h has been set to 1. 

The atoms are trapped in two internal states, treated 
here in a pseudo-spin- i framework as a = {f= \F = 
l,m F = — 1),4-= |2, — 2)}. In terms of the creation 
(a a ? ) and annihilation (a^j) operators for the respective 
spins at site i, the lattice Hamiltonian takes the two- 
component Bosc-Hubbard form: 
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where (ij) refers to nearest-neighbor sites separated by 
the lattice spacing a, U aT denotes the on-site repulsion 
between atoms in state a and r, and t a is the tunneling 
rate. The expressions for th e sp in-dependent tunneling 
energies can be found in Ref . [32| ■ For the particular case 
of 87 Rb considered here, tf « t±, Uff ~ E/4.4, ~ Uf± = U. 

The external potential V a is comprised of two terms. 
Red-detuned laser beams provide a Gaussian trap which 
is almost indistinguishable for the spin-states in consid- 
eration. In addition, there is a spin-dependent term, 
arising from a dynamically tunable external magnetic 
field, which varies linearly in space, resulting in a time- 
dependent Stern-Gerlach separation of the two spin 
states. 

The lowest excitations about the Mott insulator (U 3> 
t) do not arise from single particle tunnelings (~ t a ), but 
consist of exchanging f and 4- particles between neighbor- 
ing sites. Second order perturbation theory shows that 
the energy cost for such a process is ~ t 2 /U [32,l33|. 

In the experiment the center of the cloud is a Mott 
insulator with one particle per site. One describes this 
region via an effective Hamiltonian: 

n ef f = -jJ2s*-Sj- A v\B'\aY,iS z , (2) 

(ij) iz 

where we use spin- A operators at every site (i): Si = 

\Yjaa' a aa°aa' a a<? where Si = {S xi , S yi , S zi }. The 
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FIG. 1: (Color Online) Comparison between mean-field 
models and Monte-Carlo: (Left): magnetization (S z ) gra- 
dient as a function of /3h from Ising spin-i (solid), and clas- 
sical (dotted) mean-field theories. At \B'\ = 2G/cm, and 
T = InK, [E mag = fj,B\B'\a)/k B T ~ 15. (Right): A compar- 
ison of the energy as a function of temperature at E mag ~ J 
(\B I = O.OlG/cm), using classical mean-field theory (solid) 
and Metropolis Monte-Carlo (data points with error bars). In- 
set shows typical magnetization profiles obtained from these 
theories for the same gradient and T = InK. Position Z is 
measured in units of the lattice spacing a. The Monte-Carlo 
and mean-field curves are indistinguishable. 



Z coordinate of lattice site i is Iz- Throughout, we 
use lower case alphabets (x, y, z) to denote spin-space 
and upper case alphabets (X, Y, Z) to denote real-space. 
The super-exchange energy is J = 2t+t\/U and A/z = 
Mf ~ AH = ^Bi is the difference between the magnetic 
moments of the two states. The lattice spacing is a and 
the magnetic field gradient is B . If a different atom was 
used, one may have to correct these expressions for the 
species dependence of the tunneling, and on-site energies. 
For a lattice depth of 14.5i?R, the super-exchange scale is 
J <~ 50pK. There are two dimcnsionless numbers in the 
problem: J3J and f3E mag where E mag = A/ig\B \a and 
13 = 1/ksT. for time-independent magnetic fields, the 
Hamiltonian conserves |S | at each site, total S z and the 
total energy. 

The standard Weiss mean-field theory yields [35[ 



(S^ = t&nh(/3hi) 



(3) 



where hi = \(JY,jenn ifizj) - ApB(i)), with the sum 
over all nearest neighbor spins. Weld et al. [lj] extract 
a temperature by fitting their experimental profiles to 
Eq. |3l using the independent site approximation, J = 0. 
This technique appears to be one of the most powerful 
thermometry strategies applicable to these temperatures. 
The lowest measured temperatures in this experiment are 
T ~ 1 nK. The super-exchange energy scale is still a 
factor of 20 smaller than this, justifying the single site 
approximation. 

For this technique to provide accurate thermometry, 
one must choose the field gradient appropriately. If the 
field gradient is too large, one cannot resolve the domain 
wall. If it is too small, the finite cloud size becomes 
important. 

At the temperatures of interest, there are two contri- 
butions to the entropy density: (i) spin entropy S sp i n < 
log(2) per site and (ii) entropy from particle-hole exci- 
tations. At the final lattice depth in [lj] of 14.5Er, we 



have a t/U ~ 0.024, and a Mott gap A - 15. 5t. At 
T = InK ~ i, particle-hole excitations are suppressed 
by exp(— A/T), and at the center of the trap, nearly all 
of the entropy of the system is contained in the spin- 
excitations. Given that the Mott gap goes to zero at 
the edge of the cloud, a separate entropy reservoir exists 
here. The thermal conductivity of the Mott state is so 
low that these reservoirs can be considered decoupled. 

This poor thermal contact between the spin and parti- 
cle degrees of freedom means that the domain wall size is 
a probe of the local spin temperature and may give very 
little information of the temperature at the edge of the 
cloud. Depending on the physics one is interested in, this 
may be beneficial or detrimental. Similarly, if this "ther- 
mometer" is put in contact with another system of inter- 
est, it may be difficult for the two systems to thcrmalize 
with one another. Essentially one needs to engineer effi- 
cient mechanisms for spin exchange. Longer range spin 
interactions (such as dipolar interactions) would help in 
this regard. 



III. DYNAMICS AND THERMODYNAMICS 

Starting with the Hciscnbcrg Hamiltonian Eq. ^, we 
use the spin-commutation relations to obtain the equa- 
tions of motion for the spin operator Si = {S x , S y ,S z }. 
We then make a classical spin approximation, replacing 
these operators with their expectation value, and sequen- 
tially evolve the sites using a split-step approach. We 
verify that our timestcps are small enough that the total 
energy (for time-independent B), total S z , and |S| 2 are 
conserved to within 0.1 percent. 

In order to compare with our dynamics simulations, 
which treat the spins within a classical framework, we 
replace Eq. [3] with the the analogous mean-field theory 
for the classical model |35[ : 

M = (S zi ) = (cos(0)) = coth/% - -j-. (4) 

where hi = —2JJ2(j) $zj — A^i?(i) is the effective mag- 
netic field at site i. 

In Figfljleft) we compare the two mean-field models. 
At both small and large (3h, Eqs. [3] and 2] agree. The 
differences at intermediate f3h arise from the larger en- 
tropy of the classical model. Even though the experi- 
ments are better described by a quantum spin- ^ model, 
we work with the classical model; we do not have an 
efficient algorithm for the quantum model. At zero 
temperature, our approach can be viewed as a time- 
dependent variational ansatz, where we introduce ip = 
Ilj(cos(6j/2)exp(i(j) j /2)a\ j + sm(0 j /2)exp(-i<f>j/2)al j , 
where 9j and <fij represent the direction of the classi- 
cal spin on site j. This product ansatz should capture 
most of the gross dynamics, while drastically reducing 
the numerical complications. 

To capture thermodynamics beyond mean-field theory, 
we use a classical Metropolis Monte-Carlo algorithm. In 





FIG. 2: (Color Online) Identifying time-scales for adia- 
baticity: (Left) Time evolved magnetization profile (dots) af- 
ter Tf = Is starting with an initial field gradient of 0.25G/cm 
to a final value of O.OlG/cm at an initial temperature T = 
3nK with rate F = 20 J. The dots deviate considerably from 
the equilibrium profile from Monte-Carlo (dashed) . The solid 
lines are obtained by fitting the center and wings to two differ- 
ent temperatures using Eq(4] 0.35nK and O.lnK respectively. 
(Right) Same but with T/J — 1. This profile agrees well with 
Monte-Carlo and mean-field predictions for the equilibrium 
state at 0.35nK. Position Z is measured in units of the lattice 
spacing a. 



Fig. [TJright) we compare the classical mean-field theory 
to the Monte Carlo results. Later we compare the profiles 
obtained from the non-equilibrium dynamics to those ob- 
tained from both the classical mean-field theory and the 
Monte-Carlo. 



IV. ADIABATICITY 

Here we consider the question of exactly how slowly 
one must change parameters in order to remain in ther- 
mal equilibrium. Understanding this timescale is essen- 
tial for various proposed cooling schemes. We will focus 
on time-dependent field gradients. There are two distinct 
regimes to consider: (i) E mag ;§> J - super-exchange is 
largely irrelevant, and one requires the field gradient to 
be lowered at a rate small compared to u>l the local Lar- 
mor precession frequency and (ii) E mag < J - the gradi- 
ent should be lowered at a rate small compared to J . 

We pick a functional form for the magnetic field gradi- 
b'-b' 



ent B 



B T . where the rate T, proportional 



(l+rr) 

to B jB , characterizes the adiabaticity of the evolution: 
T/J — > oo indicates a sudden quench, while for T/J ~ 1, 
we expect the magnetization profile to adiabatically fol- 
low the equilibrium profile. 

To illustrate this we lower the magnetic field gradi- 
ent from an initial value of 0.25G/cm, (E mag ~ 25 J) to 
B \ T , = O.OlG/cm (E mag ~ J) at various rates. The ini- 
tial state is chosen to obey Eq. @] with a temperature of 
T = 3nK and with uncorrelatcd random transverse spins. 
I.E. we choose {S x ,S y } = y/1 - (S' z ) 2 {cos(0),sin(</>)}, 
where the azimuthal orientation <p is allowed to be ran- 
dom at each site. In order to more accurately reflect 
experimental conditions, we introduce spatially random 
noise in the ambient magnetic field gradient at the 10 
percent level. We ensemble average over 10 different 
noise profiles. We consider a rectangular geometry with 



Emag [in units of J] 

35 




S mag [in units of J] 

7 



|BVI [G/cm] 




0.15 
|BVI [G/cm] 



0.25 



FIG. 3: (Color Online) Cooling via adiabatic demagnetization: Left: Temperature evolution as a function of magnetic 
field gradient (\B f \(t)) and E mag = [lbB a at various times r, from \B initial \ = 2G/cm to \B F (r = 1)| = 0.25G/cm with 
T/J — 20. The temperature is extracted by fitting the spin density to the functional form in Eq. Q. Inset shows a comparison 
between the magnetization profiles from the time-evolution and Eq. [4] at T = O.lnK. Position Z is measured in units of the 
lattice spacing a. Right: Same but with \B initial \ = 0.25G/cm to \Bf inal \ = O.OlG/cm with T/J = 1, using Monte-Carlo. Lines 
through the data are guides to the eye. Inset shows the gradient of the central magnetization (oc 1/to) as a function of the 
magnetic field gradient from the time-evolution (points) and Eq. [4] (solid line) in the independent site approximation J — 0. 



5 x 5 x 30 lattice sites in the X — Y and Z— directions 
respectively. 

In Fig. [2 we plot the time-evolved magnetization pro- 
files (dots) in the direction of the field gradient, after 
ensemble averaging and averaging over the transverse di- 
rections, at Tf (~ Is) for T/J ~ 20 (left) and T/J ~ 1 
(right). By comparing the energy of the time-evolved 
profile to the classical Monte-Carlo simulation, we get 
an indication of how well thermalized the system is. For 
example, after the more rapid change in the field gra- 
dient, we find that the center of the cloud is fit better 
to a higher temperature than the wings. On the other 
hand, after the slow ramp, the coarse-grained magneti- 
zation profile and the final temperature agree well with 
the equilibrium expectations. 

Both of these ramps were slow compared to the Lar- 
mor precession frequency, so local equilibrium was main- 
tained. However the timescale for global equilibration is 
set by super-exchange. The behavior seen in Figj2] (left) 
is reminiscent of that observed by Hung et al. [36| in a 
single component Bose system, upon turning on the lat- 
tice. They too observed that different parts of the cloud 
were fit better by different temperatures. In order to 
stay close to equilibrium, it is crucial that external pa- 
rameters be changed slowly compared to the timescale 
for rearranging the state. 



ADIABATIC COOLING 



gradient adiabatically, the entropy of the system stays 
constant, and remains concentrated at the domain wall. 
This effectively means that the size of the domain wall 
does not change: the entropy contained in the domain 
wall is log(2) times the number of spins in it. Since the 
domain wall width is proportional to the ratio of the 
temperature to the magnetic energy w oc T/(Afj,^), 
one concludes that if w stays fixed, T must go down. 
Hence lowering the field gradient cools the system. Note 
that when E mag ~ J, this argument breaks down as the 
domain wall profile is set largely by the super-exchange 
energy, rather than E mag . 

In Fig. [3] (left) wc adiabatically lower the field gradient 
from 2G/cm to 0.25G/cm, at a rate of T/J = 20, start- 
ing with an initial temperature of T = InK. As we evolve 
the system, we extract the temperature by comparing the 
profiles to our Monte Carlo simulations. Due to slight de- 
viations from adiabaticity, this eight fold decrease in the 
field gradient only results in a six-fold decrease in the 
temperature. The inset shows the magnetization profile 
at the final time along with the best fit Monte-Carlo pro- 
file. 

Fig, [3] (right) illustrates what happens when E mag is 
smaller: cooling stops when E mag ~ J . The inset shows 
how the slope of the domain wall varies with field gradi- 
ent. For this simulation, we worked with a much slower 
ramp taking T = J. We found that faster ramps were 
not adiabatic. 



Having established the timescale for adiabaticity, we 
now critically evaluate how adiabatically changing the 
field gradient can lead to cooling. This is analogous to the 
adiabatic demagnetization technique used in solid state 
physics. 

Adiabatic cooling is based upon the following idea: For 
a large initial field gradient, the spin entropy is largely 
contained in the domain wall, while in the bulk the spins 
are all either pointing up or down. As we lower the field 



VI. SUPER-EXCHANGE EFFECTS AT T > J 



As shown by Quantum Monte Carlo simulations, 
l 



the 



spin- 1 Heisenberg model has a ferromagnetic transition, 
in zero field at a temperature T c ~ 0.9 J |38[ ■ Our clas- 
sical model has a similar transition at T c = 1.5J. At 
temperatures above T c , the domain wall becomes arbi- 
trarily large as B -> 0. Below T c the domain wall size 
will saturate. In addition to these gross features, there 
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FIG. 4: (Color Online) Sudden Quench: Left: Rate (■y z - 
dASz/dr) for different values of J, upon quenching from an 
initial state at T = InK, taking Bj(t) = 20(— r)G/cm. Here 
AS Z is the difference between the magnetization at the center 
and one lattice site away. Inset shows how AS Z changes in 
time, where we use St J — 0.025. Right: Typical magnetiza- 
tion profiles following the quench at r = 0, 2 J -1 , 4 J -1 , where 
flatter profiles indicate later times. Position Z is measured in 
units of the lattice spacing a. 



may be interesting physics involving the transverse spins 
in the domain wall. At high temperatures these are cer- 
tainly disordered, while at low temperature it would be 
reasonable to expect them to be ordered (at least alge- 
braically) . 

The excitations about the ferromagnet are low energy, 
quadratically dispersing, spin waves. At finite tempera- 
ture these can be thought of as non-interacting particles 
(magnons) with an average density given by the Bose oc- 
cupation factor. This picture breaks down near the Curie 
transition, as magnons can scatter off one another, and 
do not behave like free-particles. As shown by Hohenberg 
and Halperin [37| , spin- waves in the paramagnetic phase 
are over-damped, and develop a sharp peak at k = at 
T = T r . 
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FIG. 5: (Color Online) Susceptibility: Left: ASz as a func- 
tion of the field gradient at different temperatures using clas- 
sical mean- field theory: T > T c (solid), T c (dashed), and 
T « (dotted). Right: Susceptibility \ = dASz/dB '\ B , =0 in 
units of l/(G/cm) as a function of /3J from classical mean- 
field theory diverges with the mean-field exponent ~ 1.5 as 
T ^T c 



we show the spreading of the domain wall in time (r = 
0, 2J~ 1 ,4J _1 ), where the later profiles are flatter. 

The magnitude of 7 depends on the width (w) of the 
initial domain wall, which was roughly one lattice site in 
the above simulation. We find that doubling w, reduces 
7 by 2. This scaling arises because of the quadratic spin- 
wave dispersion. The group velocity v(k) = de/dk falls 
by a factor of 2 when the wave-vector is halved. We 
find that for very large initial domain walls the dynamics 
becomes extremely slow. 

These dynamics give some of the same information 
that in a traditional condensed matter setting one would 
seek from transport measurements. Such transport mea- 
surements are extremely hard in cold atoms as there is no 
natural way to introduce leads. Here one instead needs 
to infer the spin current by observing the spin dynamics. 



A. Domain wall dynamics 



B. Signatures of the phase transition 



Although the domain wall profile in Weld et al. [14[ 
does not depend on J, it is superexchange that drives 
the spin dynamics. Thus dynamics provides the clear- 
est signature of magnetic interactions in the paramag- 
netic phase. Here we explore the dynamics of the domain 
wall profile when the magnetic field gradient is suddenly 
turned off. Once the gradient is removed, the domain 
wall spreads out and flattens - presumably moving to- 
wards the equilibrium state which is a homogenous para- 
magnet. The dynamics are driven by spin currents and 
the characteristic wave-length of the relevant magnons is 
at the scale set by the domain wall size. 

We characterize the domain wall shape by the rate 
j z = dAS z /dT, where AS Z is the difference between the 
magnetization at the center, and one lattice site away. 
In Fig. H] (left) we plot this quantity as a function of J, 
starting at T = InK. We take B[{t) = 26(-T)G/cm to 
be discontinuous: 0(— r) is the Heaviside step function. 
We see that 7 oc J, this is consistent with the fact that 
J provides the scale for spin transport. In particular, 
spin wave theory would give this result. On the right 



In Fig. [5] (left) we plot AS Z as a function of the 
field gradient at different temperatures using the clas- 
sical mean-field model. The Curie transition is mani- 
fested by a divergence in the derivative of A^ with re- 
spect to the field gradient at zero field gradient (\ = 
dAS z /dB | B '_ ). On the right we plot x as a function 
of temperature, finding that as T — > T c , it diverges with 
a mean-field exponent of 1.5. 



VII. SUMMARY 

We have modeled the dynamics of the two-component 
Bose-Mott insulator by mapping it onto a spin model. 
By comparing semi-classical dynamics to classical Monte 
Carlo, and mean-field theory we extract the temperature 
of the system. We show that for experiments conducted 
in the Mott regime, the super-exchange sets the equili- 
bration rate, and adiabaticity timescales. We find that 
adiabatic demagnetization in this system can lead to sig- 
nificant cooling and expect this cooling to be observed 



experimentally. Finally the slope of the domain wall as a 
function of the field gradient is an indicator of the Curie 
transition, as it exhibits a divergence at this point. We 
propose that observing domain wall dynamics is the most 
promising avenue for studying super-exchange physics in 
the paramagnetic phase. 
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